Comprehensive identification of maize ZmE2F transcription factors and the positive role of ZmE2F6 in response to drought stress.

BACKGROUND: The early 2 factor (E2F) family is characterized as a kind of transcription factor that plays an important role in cell division, DNA damage repair, and cell size regulation. However, its stress response has not been well revealed. RESULTS: In this study, ZmE2F members were comprehensively identified in the maize genome, and 21 ZmE2F genes were identified, including eight E2F subclade members, seven DEL subfamily genes, and six DP genes. All ZmE2F proteins possessed the DNA-binding domain (DBD) characterized by conserved motif 1 with the RRIYD sequence. The ZmE2F genes were unevenly distributed on eight maize chromosomes, showed diversity in gene structure, expanded by gene duplication, and contained abundant stress-responsive elements in their promoter regions. Subsequently, the ZmE2F6 gene was cloned and functionally verified in drought response. The results showed that the ZmE2F6 protein interacted with ZmPP2C26, localized in the nucleus, and responded to drought treatment. The overexpression of ZmE2F6 enhanced drought tolerance in transgenic Arabidopsis with longer root length, higher survival rate, and biomass by upregulating stress-related gene transcription. CONCLUSIONS: This study provides novel insights into a greater understanding and functional study of the E2F family in the stress response.

Maize ZmLAZ1-3 gene negatively regulates drought tolerance in transgenic Arabidopsis.

BACKGROUND: Molecular mechanisms in response to drought stress are important for the genetic improvement of maize. In our previous study, nine ZmLAZ1 members were identified in the maize genome, but the function of ZmLAZ1 was largely unknown. RESULTS: The ZmLAZ1-3 gene was cloned from B73, and its drought-tolerant function was elucidated by expression analysis in transgenic Arabidopsis. The expression of ZmLAZ1-3 was upregulated by drought stress in different maize inbred lines. The driving activity of the ZmLAZ1-3 promoter was induced by drought stress and related to the abiotic stress-responsive elements such as MYB, MBS, and MYC. The results of subcellular localization indicated that the ZmLAZ1-3 protein localized on the plasma membrane and chloroplast. The ectopic expression of the ZmLAZ1-3 gene in Arabidopsis significantly reduced germination ratio and root length, decreased biomass, and relative water content, but increased relative electrical conductivity and malondialdehyde content under drought stress. Moreover, transcriptomics analysis showed that the differentially expressed genes between the transgenic lines and wild-type were mainly associated with response to abiotic stress and biotic stimulus, and related to pathways of hormone signal transduction, phenylpropanoid biosynthesis, mitogen-activated protein kinase signaling, and plant-pathogen interaction. CONCLUSION: The study suggests that the ZmLAZ1-3 gene is a negative regulator in regulating drought tolerance and can be used to improve maize drought tolerance via its silencing or knockout.

Research hotspots and frontiers of ethnic cultural identity--based on analysis of "web of science" database.

Cultural identity is of great significance to the formation of group consensus and the establishment of cultural self-confidence. In order to understand the history, current situation and trend, and provide theoretical support for future research, this paper makes a quantitative analysis of knowledge map including annual publication volume, trend, distribution of authors and institutions, co-occurrence, clustering and timeline of keywords as well as emergent keywords on the literature concerning ethnic cultural identity published in "Web of Science" database for a period from 2012 to 2022, with CiteSpace software as a tool. The results show an overall upward trend with diversified ethnic and regional characteristics; major institutions including universities of the U.S., the U.K., Australia, China and other countries and regions engage in their research from different disciplines such as psychology, sociology, ethnology and education; the researchers have not formed a core group of authors despite their accumulating number; research hotspots are indicated by keywords such as national identity, identity, ethnic identity and attitude; specifically, keyword clusters fall into three categories: emotional perception, multicultural identity process and ethnic cultural adaptability; researchers probe into various issues at different stages with direct relation to international situations and regional cultures. This study has positive implications for understanding and mastering the current research hotspots and development trends of ethnic cultural identity in the world.

Interkingdom multi-omics analysis reveals the effects of nitrogen application on growth and rhizosphere microbial community of Tartary buckwheat.

Tartary buckwheat (Fagopyrum tataricum Gaertn.) is an important pseudocereal crop with excellent edible, nutritional and medicinal values. However, the yield of Tartary buckwheat (TB) is very low due to old-fashioned cultivation techniques, particularly unreasonable application of nitrogen fertilizer. To improve the understanding on the theories of nitrogen use in TB, the effects of nitrogen application on growth, as well as chemical properties and microbial community of rhizosphere soil were investigated in this study. Nitrogen application could promote the plant height, stem diameter, nitrogen accumulation and yield of TB. The relative abundance and diversity of bacteria and fungi in the rhizosphere soil of TB were improved by nitrogen fertilizer. Nitrogen application increased the abundance of beneficial bacteria such as Lysobacter and Sphingomonas in rhizosphere soil, and decreased the abundance of pathogenic fungi such as Fusarium and Plectosphaerella. The results indicated that nitrogen application changed the distribution of microbial communities in TB rhizosphere soil. Furthermore, the specific enriched or depleted microorganisms in the rhizosphere soil of four TB varieties were analyzed at OTU level. 87 specific nitrogen-responsive genes with sequence variation were identified in four varieties by integrating genomic re-sequencing and transcriptome analysis, and these genes may involve in the recruitment of specific rhizosphere microorganisms in different TB varieties. This study provided new insights into the effects of nitrogen application on TB growth and rhizosphere microbial community, and improved the understanding on the mechanisms of TB root-microbe interactions.

Mapping and characterization of a novel adult-plant leaf rust resistance gene LrYang16G216 via bulked segregant analysis and conventional linkage method.

KEY MESSAGE: A novel adult-plant leaf rust resistance gene LrYang16G216 on wheat chromosome 6BL was identified and mapped to a 0.59 cM genetic interval by BSA and conventional linkage method. Leaf rust (Puccinia triticina) is one of the most devastating fungal diseases of wheat (Triticum aestivum L.). Discovery and identification of new resistance genes is essential to develop disease-resistant cultivars. An advanced breeding line Yang16G216 was previously identified to confer adult-plant resistance (APR) to leaf rust. In this research, a recombinant inbred line (RIL) population was constructed from the cross between Yang16G216 and a highly susceptible line Yang16M6393, and genotyped with exome capture sequencing and 55 K SNP array. Through bulked segregant analysis (BSA) and genetic linkage mapping, a stable APR gene, designated as LrYang16G216, was detected and mapped to the distal region of chromosome arm 6BL with a genetic interval of 2.8 cM. For further verification, another RIL population derived from the cross between Yang16G216 and a susceptible wheat variety Yangmai 29 was analyzed using the enriched markers in the target interval, and LrYang16G216 was further narrowed to a 0.59 cM genetic interval flanked by the KASP markers Ax109403980 and Ax95083494, corresponding to the physical position 712.34-713.94 Mb in the Chinese Spring reference genome, in which twenty-six disease resistance-related genes were annotated. Based on leaf rust resistance spectrum, mapping data and physical location, LrYang16G216 was identified to be a novel and effective APR gene. The LrYang16G216 with linked markers will be useful for marker-assisted selection in wheat resistance breeding.

Development of the engineered "glanded plant and glandless seed" cotton.

After fiber, cottonseed is the second most important by-product of cotton production. However, high concentrations of toxic free gossypol deposited in the glands of the cottonseed greatly hamper its effective usage as food or feed. Here, we developed a cotton line with edible cottonseed by specifically silencing the endogenous expression of GoPGF in the seeds, which led to a glandless phenotype with an ultra-low gossypol content in the seeds and nearly normal gossypol in other parts of the plants. This engineered cotton maintains normal resistance to insect pests, but the gossypol content in the seeds dropped by 98%, and thus, it can be consumed directly as food. The trait of a low gossypol content in the cottonseeds was stable and heritable, while the protein, oil content, and fiber yield or quality were nearly unchanged compared to the transgenic receptor W0. In addition, comparative transcriptome analysis showed that down-regulated genes in the ovules of the glandless cotton were enriched in terpenoid biosynthesis, indicating the underlying relationship between gland formation and gossypol biosynthesis. These results pave the way for the comprehensive utilization of cotton as a fiber, oil, and feed crop in the future.

Metabolomics and biochemical analyses revealed metabolites important for the antioxidant properties of purple glutinous rice.

Glutinous rice can be applied to many fields including brewing industry, medicine, cosmetics and food processing. However, we know very little about the basic metabolite information of glutinous rice. In this study, we identified the phenol and flavonoid metabolites in purple glutinous rice and white glutinous rice, and elucidated the relationship between metabolites and total antioxidant capacity. The results demonstrated that flavonoids contents, oligomeric proanthocyanidin contents and total antioxidant capacity of purple glutinous rice were significantly higher than those of white glutinous rice. We identified 390 differential metabolites between the purple glutinous rice and white glutinous rice by LC-MS metabolomics. Correlation analysis results showed that flavonoid and phenol metabolites contents were strongly correlated with total antioxidant capacity. This study further clarified that the pantothenic acid, pebrellin, l-glutamic acid, eupatilin, diosmin, and diosmetin could be used as candidate metabolite markers for antioxidant capacity screening in glutinous rice grains.

Zinc Transporter ZmLAZ1-4 Modulates Zinc Homeostasis on Plasma and Vacuolar Membrane in Maize.

Zinc is an essential micronutrient for plant growth and development, and functions as a cofactor for hundreds of transcription factors and enzymes in numerous biological processes. Zinc deficiency is common abiotic stress resulting in yield loss and quality deterioration of crops, but zinc excess causes toxicity for biological systems. In plants, zinc homeostasis is tightly modulated by zinc transporters and binding compounds that uptake/release, transport, localize, and store zinc, as well as their upstream regulators. Lazarus 1 (LAZ1), a member of DUF300 protein family, functions as transmembrane organic solute transporter in vertebrates. However, the function of LAZ1 in plants is still obscure. In the present study, the ZmLAZ1-4 protein was confirmed to bind to zinc ions by bioinformatic prediction and thermal shift assay. Heterologous expression of ZmLAZ1-4 in the zinc-sensitive yeast mutant, Arabidopsis, and maize significantly facilitated the accumulation of Zn2+ in transgenic lines, respectively. The result of subcellular localization exhibited that ZmLAZ1-4 was localized on the plasma and vacuolar membrane, as well as chloroplast. Moreover, the ZmLAZ1-4 gene was negatively co-expressed with ZmBES1/BZR1-11 gene through co-expression and real-time quantitative PCR analysis. The results of yeast one-hybrid and dual-luciferase assay suggested that ZmBES1/BZR1-11 could bind to ZmLAZ1-4 promoter to inhibit its transcription. All results indicated that ZmLAZ1-4 was a novel zinc transporter on plasma and vacuolar membrane, and transported zinc under negative regulation of the ZmBES1/BZR1-11 transcription factor. The study provides insights into further underlying the mechanism of ZmLAZ1-4 regulating zinc homeostasis.

Isolation and identification of a vegetative organ-specific promoter from maize.

To avoid the unregulated overexpression of the exogenous genes, specific or inducible expression is necessary for some exogenous genes in transgenic plants. But little is known about organ- or tissue-specific promoters in maize. In the present study, the expression of a maize pentatricopeptide repeat (PPR) protein encoding gene, GRMZM2G129783, was analyzed by RNA-sequencing data and confirmed by quantitative real time PCR. The results showed that the PPR GRMZM2G129783 gene specifically expressed in vegetative organs. Consequently, a 1830 bp sequence upstream of the start codon of the promoter for GRMZM2G129783 gene was isolated from maize genome (P 1830 ). To validate whether the promoter possesses the vegetative organ-specificity, the full-length and three 5'-end deletion fragments of P 1830 of different length (1387, 437, and 146 bp) were fused with glucuronidase (GUS) gene to generate promoter::GUS constructs and transformed into tobacco. The transient expression and fluorometric GUS assay in transgenic tobacco showed that all promoter could drive the expression of the GUS gene, the - 437 to - 146 bp region possessed some crucial elements for root-specific expression, and the shortest and optimal sequence to maintain transcription activity was 146 and 437 bp in length, respectively. These results indicate that the promoter of the PPR GRMZM2G129783 gene is a vegetative organ-specific promoter and will be useful in transgenic modification of commercial crops for moderate specific expression after further evaluation in monocotyledons.

Mutation of SELF-PRUNING homologs in cotton promotes short-branching plant architecture.

In cotton, the formation of fruiting branches affects both plant architecture and fiber yield. Here, we report map-based cloning of the axillary flowering mutation gene (GbAF) that causes bolls to be borne directly on the main plant stem in Gossypium barbadense, and of the clustered boll mutation gene (cl1) in G. hirsutum. Both mutant alleles were found to represent point mutations at the Cl1 locus. Therefore, we propose that the GbAF mutation be referred to as cl1b. These Cl1 loci correspond to homologs of tomato SELF-PRUNING (SP), i.e. Gossypium spp. SP (GoSP) genes. In tetraploid cottons, single monogenic mutation of either duplicate GoSP gene (one in the A and one in the D subgenome) is associated with the axillary cluster flowering phenotype, although the shoot-indeterminate state of the inflorescence is maintained. By contrast, silencing of both GoSPs leads to the termination of flowering or determinate plants. The architecture of axillary flowering cotton allows higher planting density, contributing to increased fiber yield. Taken together the results provide new insights into the underlying mechanism of branching in cotton species, and characterization of GoSP genes may promote the development of compact cultivars to increase global cotton production.

Divergence and evolution of cotton bHLH proteins from diploid to allotetraploid.

BACKGROUND: Polyploidy is considered a major driving force in genome expansion, yielding duplicated genes whose expression may be conserved or divergence as a consequence of polyploidization. RESULTS: We compared the genome sequences of tetraploid cotton (Gossypium hirsutum) and its two diploid progenitors, G. arboreum and G. raimondii, and found that the bHLH genes were conserved over the polyploidization. Oppositely, the expression of the homeolgous gene pairs was diversified. The biased homeologous proportion for bHLH family is significantly higher (64.6%) than the genome wide homeologous expression bias (40%). Compared with cacao (T. cacao), orthologous genes only accounted for a small proportion (41.7%) of whole cotton bHLHs family. The further Ks analysis indicated that bHLH genes underwent at least two distinct episodes of whole genome duplication: a recent duplication (1.0-60.0 million years ago, MYA, 0.005 < Ks < 0.312) and an old duplication (> 60.0 MYA, 0.312 < Ks < 3.0). The old duplication event might have played a key role in the expansion of the bHLH family. Both recent and old duplicated pairs (68.8%) showed a divergent expression profile, indicating specialized functions. The expression diversification of the duplicated genes suggested it might be a universal feature of the long-term evolution of cotton. CONCLUSIONS: Overview of cotton bHLH proteins indicated a conserved and divergent evolution from diploids to allotetraploid. Our results provided an excellent example for studying the long-term evolution of polyploidy.

Genomic analyses in cotton identify signatures of selection and loci associated with fiber quality and yield traits.

Upland cotton (Gossypium hirsutum) is the most important natural fiber crop in the world. The overall genetic diversity among cultivated species of cotton and the genetic changes that occurred during their improvement are poorly understood. Here we report a comprehensive genomic assessment of modern improved upland cotton based on the genome-wide resequencing of 318 landraces and modern improved cultivars or lines. We detected more associated loci for lint yield than for fiber quality, which suggests that lint yield has stronger selection signatures than other traits. We found that two ethylene-pathway-related genes were associated with increased lint yield in improved cultivars. We evaluated the population frequency of each elite allele in historically released cultivar groups and found that 54.8% of the elite genome-wide association study (GWAS) alleles detected were transferred from three founder landraces: Deltapine 15, Stoneville 2B and Uganda Mian. Our results provide a genomic basis for improving cotton cultivars and for further evolutionary analysis of polyploid crops.

Genomic insights into divergence and dual domestication of cultivated allotetraploid cottons.

BACKGROUND: Cotton has been cultivated and used to make fabrics for at least 7000 years. Two allotetraploid species of great commercial importance, Gossypium hirsutum and Gossypium barbadense, were domesticated after polyploidization and are cultivated worldwide. Although the overall genetic diversity between these two cultivated species has been studied with limited accessions, their population structure and genetic variations remain largely unknown. RESULTS: We resequence the genomes of 147 cotton accessions, including diverse wild relatives, landraces, and modern cultivars, and construct a comprehensive variation map to provide genomic insights into the divergence and dual domestication of these two important cultivated tetraploid cotton species. Phylogenetic analysis shows two divergent groups for G. hirsutum and G. barbadense, suggesting a dual domestication processes in tetraploid cottons. In spite of the strong genetic divergence, a small number of interspecific reciprocal introgression events are found between these species and the introgression pattern is significantly biased towards the gene flow from G. hirsutum into G. barbadense. We identify selective sweeps, some of which are associated with relatively highly expressed genes for fiber development and seed germination. CONCLUSIONS: We report a comprehensive analysis of the evolution and domestication history of allotetraploid cottons based on the whole genomic variation between G. hirsutum and G. barbadense and between wild accessions and modern cultivars. These results provide genomic bases for improving cotton production and for further evolution analysis of polyploid crops.

Identification and characterization of the GhHsp20 gene family in Gossypium hirsutum.

In higher plants, Heat Shock Protein 20 (Hsp20) plays crucial roles in growth, development and responses to abiotic stresses. In this study, 94 GhHsp20 genes were identified in G. hirsutum, and these genes were phylogenetically clustered into 14 subfamilies. Out of these, 73 paralogous gene pairs remained in conserved positions on segmental duplicated blocks and only 14 genes clustered into seven tandem duplication event regions. Transcriptome analysis showed that 82 GhHsp20 genes were expressed in at least one tested tissues, indicating that the GhHsp20 genes were involved in physiological and developmental processes of cotton. Further, expression profiles under abiotic stress exhibited that two-thirds of the GhHsp20 genes were responsive to heat stress, while 15 genes were induced by multiple stresses. In addition, qRT-PCR confirmed that 16 GhHsp20 genes were hot-induced, and eight genes were up-regulated under multiple abiotic stresses and stress-related phytohormone treatments. Taken together, our results presented here would be helpful in laying the foundation for understanding the complex mechanisms of GhHsp20 mediated developmental processes and abiotic stress signaling transduction pathways in cotton.

Small interfering RNAs from bidirectional transcripts of GhMML3_A12 regulate cotton fiber development.

Natural antisense transcripts (NATs) are commonly observed in eukaryotic genomes, but only a limited number of such genes have been identified as being involved in gene regulation in plants. In this research, we investigated the function of small RNA derived from a NAT in fiber cell development. Using a map-based cloning strategy for the first time in tetraploid cotton, we cloned a naked seed mutant gene (N1 ) encoding a MYBMIXTA-like transcription factor 3 (MML3)/GhMYB25-like in chromosome A12, GhMML3_A12, that is associated with fuzz fiber development. The extremely low expression of GhMML3_A12 in N1 is associated with NAT production, driven by its 3' antisense promoter, as indicated by the promoter-driven histochemical staining assay. In addition, small RNA deep sequencing analysis suggested that the bidirectional transcriptions of GhMML3_A12 form double-stranded RNAs and generate 21-22 nt small RNAs. Therefore, in a fiber-specific manner, small RNA derived from the GhMML3_A12 locus can mediate GhMML3_A12 mRNA self-cleavage and result in the production of naked seeds followed by lint fiber inhibition in N1 plants. The present research reports the first observation of gene-mediated NATs and siRNA directly controlling fiber development in cotton.

Genetic basis for glandular trichome formation in cotton.

Trichomes originate from epidermal cells and can be classified as either glandular or non-glandular. Gossypium species are characterized by the presence of small and darkly pigmented lysigenous glands that contain large amounts of gossypol. Here, using a dominant glandless mutant, we characterize GoPGF, which encodes a basic helix-loop-helix domain-containing transcription factor, that we propose is a positive regulator of gland formation. Silencing GoPGF leads to a completely glandless phenotype. A single nucleotide insertion in GoPGF, introducing a premature stop codon is found in the duplicate recessive glandless mutant (gl2gl3). The characterization of GoPGF helps to unravel the regulatory network of glandular structure biogenesis, and has implications for understanding the production of secondary metabolites in glands. It also provides a potential molecular basis to generate glandless seed and glanded cotton to not only supply fibre and oil but also provide a source of protein for human consumption.

Sequencing of allotetraploid cotton (Gossypium hirsutum L. acc. TM-1) provides a resource for fiber improvement.

Upland cotton is a model for polyploid crop domestication and transgenic improvement. Here we sequenced the allotetraploid Gossypium hirsutum L. acc. TM-1 genome by integrating whole-genome shotgun reads, bacterial artificial chromosome (BAC)-end sequences and genotype-by-sequencing genetic maps. We assembled and annotated 32,032 A-subgenome genes and 34,402 D-subgenome genes. Structural rearrangements, gene loss, disrupted genes and sequence divergence were more common in the A subgenome than in the D subgenome, suggesting asymmetric evolution. However, no genome-wide expression dominance was found between the subgenomes. Genomic signatures of selection and domestication are associated with positively selected genes (PSGs) for fiber improvement in the A subgenome and for stress tolerance in the D subgenome. This draft genome sequence provides a resource for engineering superior cotton lines.

The R3-MYB gene GhCPC negatively regulates cotton fiber elongation.

Cotton (Gossypium spp.) fibers are single-cell trichomes that arise from the outer epidermal layer of seed coat. Here, we isolated a R3-MYB gene GhCPC, identified by cDNA microarray analysis. The only conserved R3 motif and different expression between TM-1 and fuzzless-lintless mutants suggested that it might be a negative regulator in fiber development. Transgenic evidence showed that GhCPC overexpression not only delayed fiber initiation but also led to significant decreases in fiber length. Interestingly, Yeast two-hybrid analysis revealed an interaction complex, in which GhCPC and GhTTG1/4 separately interacted with GhMYC1. In transgenic plants, Q-PCR analysis showed that GhHOX3 (GL2) and GhRDL1 were significantly down regulated in -1-5 DPA ovules and fibers. In addition, Yeast one-hybrid analysis demonstrated that GhMYC1 could bind to the E-box cis-elements and the promoter of GhHOX3. These results suggested that GhHOX3 (GL2) might be downstream gene of the regulatory complex. Also, overexpression of GhCPC in tobacco led to differential loss of pigmentation. Taken together, the results suggested that GhCPC might negatively regulate cotton fiber initiation and early elongation by a potential CPC-MYC1-TTG1/4 complex. Although the fibers were shorter in transgenic cotton lines than in the wild type, no significant difference was detected in stem or leaf trichomes, even in cotton mutants (five naked seed or fuzzless), suggesting that fiber and trichome development might be regulated by two sets of genes sharing a similar model.

Anthropogenic halo disturbances alter landscape and plant richness: a ripple effect.

Although anthropogenic landscape fragmentation is often considered as the primary threat to biodiversity, other factors such as immediate human disturbances may also simultaneously threaten species persistence in various ways. In this paper, we introduce a conceptual framework applied to recreation landscapes (RLs), with an aim to provide insight into the composite influences of landscape alteration accompanying immediate human disturbances on plant richness dynamics. These impacts largely occur at patch-edges. They can not only alter patch-edge structure and environment, but also permeate into surrounding natural matrices/patches affecting species persistence-here we term these "Halo disturbance effects" (HDEs). We categorized species into groups based on seed or pollen dispersal mode (animal- vs. wind-dispersed) as they can be associated with species richness dynamics. We evaluated the richness of the two groups and total species in our experimental landscapes by considering the distance from patch-edge, the size of RLs and the intensity of human use over a six-year period. Our results show that animal-dispersed species decreased considerably, whereas wind-dispersed species increased while their richness presented diverse dynamics at different distances from patch-edges. Our findings clearly demonstrate that anthropogenic HDEs produce ripple effects on plant, providing an experimental interpretation for the diverse responses of species to anthropogenic disturbances. This study highlights the importance of incorporating these composite threats into conservation and management strategies.

Inheritance of long staple fiber quality traits of Gossypium barbadense in G. hirsutum background using CSILs.

Gossypium hirsutum is a high yield cotton species that exhibits only moderate performance in fiber qualities. A promising but challenging approach to improving its phenotypes is interspecific introgression, the transfer of valuable traits or genes from the germplasm of another species such as G. barbadense, an important cultivated extra long staple cotton species. One set of chromosome segment introgression lines (CSILs) was developed, where TM-1, the genetic standard in G. hirsutum, was used as the recipient parent and the long staple cotton G. barbadense Hai7124 was used as the donor parent by molecular marker-assisted selection (MAS) in BC(5)S(1­4) and BC(4)S(1­3) generations. After four rounds of MAS, the CSIL population was comprised of 174 lines containing 298 introgressed segments, of which 86 (49.4%) lines had single introgressed segments. The total introgressed segment length covered 2,948.7 cM with an average length of 16.7 cM and represented 83.3% of tetraploid cotton genome. The CSILs were highly varied in major fiber qualities. By integrated analysis of data collected in four environments, a total of 43 additive quantitative trait loci (QTL) and six epistatic QTL associated with fiber qualities were detected by QTL IciMapping 3.0 and multi-QTL joint analysis. Six stable QTL were detected in various environments. The CSILs developed and the analyses presented here will enhance the understanding of the genetics of fiber qualities in long staple G. barbadense and facilitate further molecular breeding to improve fiber quality in Upland cotton.